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SUMMARY 


A new incompressible lifting-surface theory is developed for thin rectan- 
gular wings. The solution requires the downwash equation to be in the form of 
Cauchy-type integrals. Lan T s method is employed for the chordwise integrals 
since it properly accounts for the leading-edge singularity, Cauchy singularity 
and Kutta condition. The Cauchy singularity in the spanwise integral is also 
accounted for by using the midpoint trapezoidal rule and theory of Chebychev 
polynomials. The resulting matrix equation, formed by satisfying the boundary 
condition at control points, is simplier and quicker to compute than other lift- 
ing surface theories. Solutions were found to converge with only a small num- 
ber of control points and to compare favorably with results from other methods. 


INTRODUCTION 


Numerous subsonic lifting surface theories have been developed for thin 
wings over the past thirty years. A comparison of three of the more prominent 
numerical methods in 1968 is given in reference 1. Although all these methods 
give essentially the same results, they require tedious integration techniques 
and consume considerable computational time. Vortex lattice methods are sim- 
pler and can be applied to more complex configurations; however, they are gen- 
erally less accurate than lifting surface methods. This author developed a 
lifting surface method for thin rectangular wings (ref. 2) which can also be 
interpreted as a vortex lattice method. The present paper examines the conver- 
gence of solutions and compares results with those of the NRL (National Aero- 
space Laboratory, Netherlands) method given in reference 1. 


SYMBOLS 


A aspect ratio 

b wing span 

c wing chord 


* This research is supported by the U. S. Army Research Office, Research Tri- 
angle Park, N. C., under Grant Number DAAG29-76-G-0045 . 
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sectional lift coefficient 
sectional drag coefficient 
far-field induced-drag coefficient 

near-field induced - * drag coefficient 

wing lift coefficient 

wing pitching-moment coefficient, about leading edge 

leading-edge suction parameter, see eq. (13) 

parameter defined by eq. (10) 
parameter defined by eq. (15) 

number of spanwise control points over whole span 
number of chordwise control points 
National Aerospace Laboratory, Netherlands 
wing planform area 
freestream velocity 


w 


x 


non-dimensional downwash velocity, referred to and positive up- 
wards 

chordwise coordinate measured from leading edge in direction of V 

00 


ac 


X 


ac 


sectional and wing aerodynamic-center locations, respectively 


y 

a 

Y 

r 

0 


spanwise coordinate, positive to the right 
angle of attack 

non-dimensional circulation per unit chord 
circulation 

transformed chordwise coordinate, see eq. (7) 
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(j> transformed spanwise coordinate, see eq. (8) 

Subscripts: 


i 

3 

k 

a 

p 


chordwise control point, see eq. (5) 
spanwise control point, see eq. (17) 
chordwise integration point, see eq. (4) 
spanwise integration point, see eq. (16) 
evaluated at spanwise position (j) = pTt/M 


ANALYSIS 


For simplicity, the present method is developed for rectangular 
wings. The downwash equation from lifting surface theory is usually given by 
one of the two following forms (ref. 3): 
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Equation (1) is the form used by the three methods compared in reference 1. It 
contains the Mangier-type integral due to the term (y - y^) 2 in the denominator. 
Equation (2), however, contains Cauchy-type integrals since the terms (y - y^) 
and (x - X]J in the denominator are linear. The present method requires the 
downwash to be in the form of equation (2) . 


In order to understand the development of the present method, con- 
sider the two-dimensional problem first. Lan (ref. 4) developed an ingenious 
method for thin airfoils by using the midpoint trapezoidal rule and the theory 
of Chebychev polynomials to reduce the two-dimensional downwash integral 
(Cauchy- type integral) to a finite sum. The remarkable accuracy of this tech- 
nique is due to the following summational result: 


N 


l 

k=l 


1 


COS 0- - cos 0 . 
k x 


r-r 


L N 


( i = o) 

( i + 0,N) 

( i = n) 
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where 


(4) 


n _ (2k - 1 )tt 
0 k 2N 


(k 1, • • • » N ^ 


are the integration points and 

0. ® i7r/N 
i 



(5) 


are the control points (positions where the boundary condition is applied) . 
Note the similarity of equation (3) to the integral result 
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cos 0^ - cos 0 


= 0 


( 6 ) 


Lan used these equations to develop the two-dimensional downwash summation which 
properly accounts for the Cauchy singularity, the leading-edge square-root sin- 
gularity, and the Kutta condition at the trailing edge. 

Since equation (2) contains Cauchy-type integrals, Lan's approach 
will be applied to both the chordwise and spanwise integrals. First, transform 
the chordwise coordinate by 

2 x/c = 1 - cos 0 (7) 

and the spanwise coordinate by 

2y/b = - cos (J) (8) 

Then, use Multhopp T s interpolation formula to represent the circulation per 
unit chord by the trigonometric sum 

M-l M-l 

Y(0»4>) = I Y (9) I sin neb sin n<J> (9) 

p=l F n=l p 

where the subscript "p" refers to the spanwise position (j>p = pm/M and (M - 1) 
represents the number of spanwise control points. Substitute equations (7), 

(8), and (9) into equation (2) and consider the chordwise integral first. With 

G(0,0^,cj),(j)^) E /(cos 0^ - cos + A^(sin <J>^ - sin (f))^ (10) 

the chordwise integrals are 
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The integral on the right side can be reduced to a finite sum by the midpoint 
trapezoidal rule; and equations (3) and (6) can be used to account for the 
leading-edge singularity, Cauchy singularity, and Kutta condition (y p (tt) = 0) as 
shown below: 


J Y p^ 6 l^ sin 9 1 G d0 l 
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where 0^. are the chordwise integration points given by equation (4) and 0^ are 
the chordwise control points given by equation (5) . The leading-edge suction 
parameter is defined as 


4C„ = lim y (0) sin 0 

S p 9+0 p 


(13) 


Now substitute equation (12) into equation (2) and perform the spanwise inte- 
gration in a somewhat similar manner, accounting for the Cauchy singularity, to 
obtain the final form of the downwash as 
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The spanwise integration points are 

, (2J£ - 1) TT 

h 2M 


(s, 1, • • • » M ) 
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and the spanwise control points are 


4>j = jtt/m 


( j = 1 M - 1 ) 

The tangent-flow boundary condition for flat wings requires that w i j = 


(17) 


- a. 
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The N(M - 1) values of Yp s k are calculated by solving the matrix equation formed 
by applying equation (14; for i ^ 0 at the chordwise and spanwise control points 
given by equations (5) and (17). Then after the Yp 5 k are calculated, the (M -1) 
leading-edge suction parameters Cgj can be computed by successively applying 
equation (14) with i - 0 (control point at the leading edge) at the spanwise 
positions j = 1, . .., M - 1. Regardless of the number (N) of chordwise control 
points used, there is always a control point at the trailing edge which accounts 
for the Kutta condition, and another control point at the leading edge which 
gives the leading-edge suction parameter, if desired. 

The sectional and wing aerodynamic characteristics may now be cal- 
culated by using the midpoint trapezoidal rule to reduce the integrals to fi- 
nite sums, as illustrated below: 
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The spanwise loading can be made continuous by equation (9) , and the chordwise 
loading can also be made continuous by fitting Cgj and Yj ^ to the chordwise 
loading functions for thin airfoil theory. 


RESULTS AM) DISCUSSION 


For one control point (N = 1, M - 2) the present method yields 

CL /a = ttA/ (1 + A + A 2 /2) and CL = C^/ttA 
L D . L 

1 

These results give the correct limit as A -*■ 0, but just as Lan found for air- 
foils, at least two chordwise control points are needed to get an accurate 
pitching moment. Reference 2 showed that the present spanwise integration 
method gives the exact classical solution to Prandtl f s lifting line equation. 
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Table 1 gives a detailed comparison of the present method with the 

NLR method (ref. 1) for an A = 2 rectangular wing with N = 4 and M = 16. Al- 

though the NLR method used 15 spanwise loading functions, 127 spanwise integra- 
tion points were employed. Excellent agreement is obtained between the two 
methods for the lift, pitching moment, aerodynamic center, far-field drag, and 
spanwise lift distribution. The spanwise variation of section drag and aero- 
dynamic center compares well except near the wing tips, and the near-field drag 
values differ. As noted in reference 1, at least 8 chordwise control points 
are needed to get convergence of the section drag near the wing tips and the 

near-field drag. The computational time required for the results in Table 1 

was 22 minutes for the NLR method on a CDC 3300 computer, whereas the present 
method required less than 10 seconds on an IBM 370/165 computer. 

The effects of the number of control points on the convergence of 
lift and aerodynamic center are shown in figure 1 for rectangular wings with 
A = 2 and 7. These results indicate that good results are obtained for the 
A = 7 wing with N = 2 and M = 10, whereas more chordwise and less spanwise con- 
trol points are needed for the A = 2 wing. Figure 2 illustrates the conver- 
gence of the near- and far-field induced drag for the same two wings. Note 
that the far-field induced drag is insensitive to both M and N for the A = 2 
and A = 7 wings. On the other hand, the near-field induced drag depends on both 
M and N, particularly for the A = 2 wing. 


CONCLUDING REMARKS 


The present lifting-surface method for rectangular wings was found 
to compare favorably with other methods, but it is simpler and requires smaller 
computational times. The number of control points required for convergence of 
the aerodynamic characteristics is dependent on both the wing aspect ratio and 
the aerodynamic parameter. Convergence is fast for lift, pitching moment, 
aerodynamic center, far-field drag, and spanwise lift distribution. For the 
A = 7 wing two chordwise and about 10 spanwise control points gave good results, 
whereas at least 4 chordwise and 8 spanwise control points are needed for the 
A = 2 wing. The far-field induced drag is particularly insensitive to the num- 
ber of control points, with good results for N = 2 and M = 4 on both wings. 
Convergence is slow for the section drag and aerodynamic center near the wing 
tips. The near-field drag converges more slowly than any of the other para- 
meters. Other planforms are presently being investigated. 
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TABLE 1. RESULTS FOR RECTANGULAR PLANFORM, A = 2 


Overall Values 



Present 

NLR 

(ref. 1) 


N=4, M=16 

N=4 , M=16 

C L /a 

2.4732 

2.4744 

-v® 

0.5187 

0.5182 

X /c 
ac 

0.2097 

0.2094 

” ac d. /c l 

1 

1.0007 

1.0007 

11 

0.9951 

1.0108 


Values of 



2y/b 

Present 

NLR 

0 

0.1847 

0.1848 

0.1951 

0.1832 

0.1832 

0.3827 

0.1784 

0.1781 

0.5556 

0.1693 

0.1686 

0.7071 

0.1548 

0.1541 

0.8315 

0.1331 

0.1353 

0.9239 

0.0988 

0.1131 

j 0.9808 

0.0394 

0.0770 


Values of c /C 

Li 


Values of x /c 
ac 


2y/b 

Present 

NLR 


2y/b 

Present 

NLR 

0 

1.2543 

1.2543 


0 

0.2200 

0.2199 

0.1951 

1.2331 

1.2331 

• 

0.1951 

0.2187 

0.2187 

0.3827 

1.1692 

1.1692 


0.3827 

0.2150 

0.2149 

0.5556 

1.0625 

1.0625 


0.5556 

0.2087 

0.2085 

0.7071 

0.9137 

0.9137 


0.7071 

0.1999 

0.1996 

0.8315 

0.7257 

0.7257 


0.8315 

0.1896 

0.1886 

0.9239 

0.5045 

0.5044 


0.9239 

0.1798 

0.1773 

0.9808 

0.2588 

0.2587 


0.9808 

0.1731 

0.1685 
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